Excitonic effects in the optical properties of SiC sheet and nanotubes 
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The quasiparticle band structure and optical properties of single-walled zigzag and armchair SiC 
nanotubes (SiC-NTs) as well as single SiC sheet are investigated by ab initio many-body calculations 
using the GW and the GW plus Bethe-Salpeter equation (GW+BSE) approaches, respectively. 
Significant GW quasiparticle corrections of more than 1.0 eV to the Kohn-Sham band gaps from 
the local density approximation (LDA) calculations are found. The GW self-energy corrections 
transform the SiC sheet from a indirect LDA band gap to a direct band gap material. Furthermore, 
the quasiparticle band gaps of SiC-NTs with different chiralities behave very differently as a function 
of tube diameter, and this can be attributed to the difference in the curvature- induced orbital 
rehybridization between the different chiral nanotubes. The calculated optical absorption spectra 
are dominated by discrete exciton peaks due to exciton states with large binding energy up to 
2.0 eV in the SiC sheet and SiC-NTs. The formation of strongly bound excitons is attributed to 
the enhanced electron-hole interaction in these low dimensional systems. Remarkably, the excited 
electron amplitude of the exciton wavefunction is found to peak on the Si atoms near the hole position 
(which is on the C site) in the zigzag SiC-NTs, indicating a charge transfer from an anion (hole) to 
its neighboring cations by photoexcitation. In contrast, this pronounced peak structure disappear 
in the exciton wavefunction in the armchair SiC-NTs. Furthermore, in the armchair SiC-NTs, the 
bound exciton wavefunctions are more localized and also strongly cylindrically asymmetric. The 
large excitation energy of ^ 3.0 eV of the first bright exciton with no dark exciton below it, suggests 
that the small-radius armchair SiC-NTs be useful for optical devices working in the UV regime. 
On the other hand, the zigzag SiC-NTs have many dark excitons below the first bright exciton and 
hence may have potential applications in tunable optoelectric devices ranging from infrared to UV 
frequencies by external perturbations. 

PACS numbers: 81.07.De, 73.22.-f, 78.67.Ch,73.63.Fg 



INTRODUCTION 

Silicon carbide (SiC) crystallizes in either a cubic 
or a hexagonal form, and exbihits interesting polytyp- 
ism [l|, m . The polytypes are made of identical hexagonal 
layers with different stacking sequences. These polytypes 
are semiconductors with a range of band gaps, from 2.39 
eV in the zincblende polytype (3C) to 3.33 eV in the 
wurtzite polytype (2H) [l|, l3|. Furthermore, 3C and 6H 
SiC are used for hi^hjbemperature, high-power and hi^- 
frequency devices [4- 7] due to their unique properties [8] , 
while 6H SiC with a band gap of 2.86 eV is a useful ma- 
terial for blue light-emitting diode applications [3]. 

Recently, SiC nanotubes (NTs) were synthesized by re- 
action of carbon nanotubes (CNTs) with SiO at various 
temperatures p. Like other tubular materials such as 
C [iMi, BN [l3,[i3, AIN [Isl, and GaN nanotubes 
that have been synthesized previously, SiC-NTs display 
very interesting properties distinctly different from their 
bulks. Therefore, the successful synthesis of the SiC-NTs 
has stimulated a number of theoretical and experimen- 
tal investigations on tubular form of SiC (see [13, [lH 
and references therein). In particular, based on density- 
functional calculations, Miyamoto and Yu [19] predicted 



that the strain energies of SiC-NTs are lower than that 
of CNTs, and that the band gaps of SiC-NTs can be di- 
rect or indirect, depending on the chirality. Using both 
tight-binding molecular dynamics and ab initio methods, 
Menon and coworkers [1^ showed that single- walled SiC- 
NTs are very stable with a large band gap. 

Unlike CNTs, SiC-NTs are polar materials and there- 
fore, may exhibit unusual physical properties that CNTs 
may not have. For example, like BN-NTs (see, e.g., 
0,|22|), zigzag SiC-NTs may become piezoelectric, and 
also show second-order non-linear optical response. A 
knowledge of the optical properties of SiC-NTs is im- 
portant for their optical and electro-optical applications. 
Therefore, Wu and Guo have recently carried out a series 
of ab initio calculations within the density functional the- 
ory (DFT) with the local density approximation (LDA) 
in order to analyze the linear optical features and under- 
lying band structure as well as the second-harmonic gen- 
eration and linear electro-optical coefficients of all three 
types of the SiC-NTs. [l7|, [isl In particular, Wu and Guo 
found that all the SiC nanotubes are semiconductors with 
exceptions of the ultrasmall (3,0) and (4,0) zigzag tubes 
which may be regarded as the thinnest conducting SiC 
nanowires.[l7j Interestingly, the band gap of the zigzag 
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SiC-NTs may be reduced from the energy gap of the SiC 
sheet ah the way down to zero by reducing the diameter, 
though the band gap for ah the SiC nanotubes with a 
diameter larger than ~20 A approaches that of the SiC 
sheet. Furthermore, DFT indicates that all the semi- 
conducting zigzag SiC-NTs have a direct band gap. All 
these suggest that they may have interesting applications 
in optical and optoelectronic devices. Nonetheless, both 
the armchair and chiral SiC-NTs have an indirect band 
gap at the DFT level. 

These previous theoretical studies of the electronic 
and optical properties of SiC-NTs were based on the 
independent-particle approximation within the DFT 
framework [itI Hsj or on semi-empirical calculations |2q|. 
However, many-body interactions are known to play 
an important role in low-dimensional systems such as 
nanotubes I23L 124]. due to reduced charge screening and 
enhanced electron-electron correlation. Therefore, moti- 
vated by the prospects of optoelectronic device applica- 
tions of SiC-NTs and also this theoretical issue, we have 
performed ab initio calculations to study the many-body 
effects on the quasiparticle band gaps and optical spectra 
of SiC-NTs. Indeed, we find pronounced GW quasipar- 
ticle corrections of more than 1.0 eV to the LDA band 
gaps in the SiC-NTs. Further, the calculated optical ab- 
sorption spectra from the SiC-NTs are dominated by the 
discrete exciton peaks due to strongly bound excitons 
with a large binding energy of up to 2.0 eV. 

SiC-NTs can be considered as an atomic layer of SiC 
sheet rolled up into a cylinder, and similar to the CNTs 
the structure of a SiC-NT is specified by the chiral vector 
which is given in term of a pair of integers (n, m).[ll] Un- 
derstanding the optical properties of an isolated SiC sheet 
would help understand the many-body effects on the op- 
tical responses of SiC-NTs. Therefore, the single atomic 
layer hexagonal SiC sheet is first considered here. Fur- 
thermore, the optical properties of the two-dimensional 
(2D) single SiC sheet are interesting on their own ac- 
count, and, in particular, a single SiC sheet bears many 
similarities to graphene which exhibits many fascinating 
properties. [25, 26] 

This paper is organized as follows. In the next sec- 
tion, the theoretical methods and computational details, 
including the ground-state, quasiparticle and electron- 
hole interaction calculations, are described. In Sec. Ill, 
the calculated quasipartice band structure and optical 
absorption spectrum of the isolated SiC sheet are pre- 
sented. In Sec. IV, the calculated quasiparticle band 
structures of the SiC-NTs are reported and the effects of 
curvature and chirality on the band gaps are analyzed. 
In Sec. V, the calculated optical absorption spectra of 
the SiC-NTs are presented, and the excitonic effects are 
discussed. Finally, the conclusions drawn from this work 
are given in Sec. VI. 



METHODOLOGY 

In the present paper, we adopt the approach of Rohlf- 
ing and Louie [271, l28| and hence calculate the quasi- 
particle energies, electron-hole excitations and optical 
spectra from first-principles in three steps. First, we 
calculate the electronic ground state with DFT in the 
LDA.fii^ Second, we obtain the quasiparticle energies 
(E^^) within the GW approximation for the electron 
self energy i;.[^, (sij Finally, we evaluate the coupled 
electron-hole excitation energies and optical spectra by 
solving the Bethe-Salpeter equation (BSE) of the two- 
particle Green's function. |27l. l28| Details of our compu- 
tations are described below. 



Ground-state properties 

An appropriate mean-field description of the ground- 
state properties of solids is essential to perform quasipar- 
ticle calculations within many-body perturbation theory. 
For the conventional covalent and ionic bonded materials, 
it is shown |3l| that eigen wavefunctions and eigen ener- 
gies of the Kohn-Sham equation [29] in the LDA provide a 
good starting point for the many-body perturbation cal- 
culations such as within the GW approximations. There- 
fore, in this paper, the ground-state electronic struc- 
ture of the single SiC sheet and SiC-NTs are first calcu- 
lated by means of ab initio plane-wave pseudopotential 
methodfsij. The electronic configurations of Si(35^,3p^) 
and C(2s^, 2p'^) are treated as active valence states when 
generating pseudopotentials in the Kleinman-Bylander 
form [33] for all the calculations. The cutoff energy of 40 
Ry for the plane wave basis set is used. Also to ensure 
the convergency, the Brillouin zone (BZ) integrals for the 
calculations of the SiC sheet and SiC-NTs are carried 
out using the 10x10x1 and 1x1x32 Monkhorst-Pack k- 
grid sampling [34], respectively. The supercell approach 
is used here such that an isolated SiC sheet (nanotube) 
is approximated by a single SiC layer (nanotube) sur- 
rounded by vacuum in a supercell. Sufficient intersheet 
and intertube separations (more than 10 A) are used to 
prevent intersheet and intertube interactions. The un- 
derlying atomic positions and lattice constants are taken 
from Ref. fl^ where the theoretical structures were ob- 
tained by the conjugate gradient method with the atomic 
forces and the stress calculated from first-principles. 



Quasiparticle calculations 

Based on the ground-state Kohn-Sham wavefunctions 
and corresponding eigenvalues (E^^) calculated above, 
the many-body effects on the quasiparticle band struc- 
ture characterized by the self energy (S) can be evaluated 
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FIG. 1. The calculated quasiparticle band gap of a SiC sheet 
(open symbols) as a function of the inverse of intersheet sep- 
aration (Lz)- The converged band gap indicated as the solid 
line is derived by employing the Coulomb truncation scheme. 
The dashed line linking the open circles is a guide to the eye 
only. 



by solving the Dyson equation 



Ho + S(E«f ) 



\n\^>=E^^\nk> 



(1) 



where Hq is the Hamiltonian in the Hatree approxima- 
tion, and |nk > represents the quasiparticle wavefunc- 
tion of energy within one-partcile Green's function 
method. In the G^W^ approximation, the vertex cor- 
rection is approximated by a delta function and the irre- 
ducible polarizibility is a convolution of the mean field 
Green's function G^, i.e., = —iG^G^. This gives rise 
to a dynamical dielectric matrix within the random phase 
approximation (RPA) scheme and a generalized plasmon- 
pole (GPP) model[3l|. Also, the screened Coulomb po- 
tential can be obtained from the bare Coulomb in- 
teraction V through the equation = v / \l — vP^'\. 
Finally, the quasiparticle self-energy operator is given by 
E = iG^W^. Although alternatives have been proposed 
recently to calculate the self energy going beyond this 
approximation, the G'^W'^ approach is still the most ef- 
ficient scheme for providing a successful description of 
the electronic excitation and even transport properties 
in various semiconductors, metals, surfaces, and novel 
nanomaterials [36l-[38| . Therefore, in this paper, we adopt 
the G^W'^ approach to calculate the quasiparticle prop- 
erties of the single SiC sheet and SiC-NTs. Here a rather 
dense /c-point sampling for the SiC sheet (18x18x1) and 
also for the SiC-NTs (1x1x32) is used to ensure the cal- 
culated quasiparticle energies converged to within 0.05 
eV. 

A proper description of the dynamical screening prop- 



erties of a solid plays an important role in obtaining reli- 
ably its quasiparticle properties. In low dimensional sys- 
tems, such requirement becomes crucial due to the highly 
anisotropic screening interactions associated mainly with 
the vacuum region between the isolated units in the su- 
per cell scheme. A larger vacuum region is preferred to 
reduce the residual unphysical inter-unit interaction and 
eventually ends up with an ideal vacuum screening of 
the isolated units. For example, the band gap of an 
isolated SiC sheet, as shown in Fig. 1, would increase 
as the screening from neighboring layers is reduced by 
increasing the separation {Lz) between the SiC layers 
in the neighboring supercells. Furthermore, the conver- 
gence of the band gap with respect to the spatial sepa- 
ration of the neighboring SiC layers is very slow. Due 
to the computation constraints, it is impossible to ob- 
tain the fully converged band gap simply by expanding 
the inter-layer separation. Therefore, an efficient trunca- 
tion scheme [H, liol to the Coulomb interaction in which 
a step function at the boundary of the supercell is in- 
troduced, will be adopted in the present GW and BSE 
calculations to remove the long-range Coulomb interac- 
tion between the structure and its images. Fig. 1 clearly 
demonstrates that the truncation scheme indeed offers 
an efficient means to achieve the converged quasiparti- 
cle band gap for the low-dimensional structures. Table 
I further shows that the quasiparticle band gap and also 
quasiparticle energy differences from the GW calcula- 
tions with or without the Coulomb truncation scheme 
can differ significantly. 



Electron-hole excitations 

The optical properties of a sohd, associated with the 
interaction between light and the electronic excitations 
of a system, are described by the frequency-dependent 
macroscopic dielectric function e{uj). For example, the 
optical absorption spectrum is determined by the imagi- 
nary part of the dielectric function, 62 (cj). According to 



TABLE I. The band gap (Egap) and quasiparticle energy dif- 
ferences between the conduction band minimum (c) and va- 
lence band maximum (v) at some high symmetry /c-points 
of the SiC sheet from the LDA and GW [without/with the 
Coulomb truncation (no trunc.)/(trunc.)] calculations. All 
energies are in eV. 
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4.46 
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the Fermi Golden rule, the transition energy associated 
with a peak in the optical absorption spectrum can be 
estimated from the energy difference between the associ- 
ated electron and hole states. The corresponding oscil- 
lator strength can be obtained from the optical transi- 
tion matrix elements derived from the electronic ground 
and excited states involved. Since the incident photons 
carry negligible momentum, the peaks in the optical ab- 
sorption spectrum can be considered as direct interband 
transitions between the occupied and unoccupied elec- 
tronic states without any crystal momentum transfer. 
Therefore, if ignoring the electron-hole interaction, both 
linear and nonlinear optical spectra can be calculated 
in an independent-particle approximation by including 
only the vertical interband transitions between the Kohn- 
Sham states [13, [HI. In the GW+RPA approach, the 
dielectric function calculated based on the quasiparticle 
energies as described above, gives the interband optical 
absorption spectrum beyond the LDA, but still without 
electron-hole interaction (or excitonic) effects. 

However, for the systems with strong excitonic effects, 
the optical properties may be dominated by excitons 
which are composed of strongly correlated electron-hole 
pairs of the systems. The connection of the exciton ener- 
gies and corresponding electron-hole amplitudes A^^^ 
of the correlated electron-hole excitations S is governed 
by the Bethe-Salpeter equation[28] 



(2) 

where ^vk(^ck) denotes the quasiparticle (e.g., GW) 
eigenvalues of valence (conduction) bands at a specific 
k-point and K^^ is the kernel describing the interaction 
between excited electrons and holes. Therefore, a more 
realistic two-particle picture of the optical excitations in- 
cluding excitonic effects in these systems should be pro- 
vided by means of the GW+BSE method. To obtain 
the converged optical spectra, in the present work, the 
kernel K""^ of the SiC sheet (SiC-NTs) is computed on 
a sparse /c-point grid of 18x18x1 (1x1x32) and then 
interpolated [l^ onto a denser /c-point grid of 36x36x1 
(1x1x64). 

Instructively, the electron-hole amplitude in real space 
can be expanded in the quasielectron and quasihole basis 
{0ck(re),^^k(rh)}, 



^s{re^rh) = ^^A^ck</>ck(re)(/>^k(r^) 



(3) 



k c,v 



and the corresponding exciton states can be visualized 
in real space. Because of the complexity of the six- 
dimensional exciton wavefunction, a simpler distribu- 
tion of the electron amplitude square with the hole po- 
sition fixed, i.e., \^s{'^e^'^h = 0)1^, is usually used to 
reveal the essence of the electron-hole correlation of the 



exciton. [23|,|2J 

Finally, as noticed before 12. [l7| . the employment of 
a 3D supercell method for reduced-dimensional systems 
will generate an arbitrary volume effect on the dielec- 
tric function computation. In order to resolve this am- 
biguity and compare directly to experimental measure- 
ments, we calculate the imaginary part of polarizability 
per unit area (in units of nm) for the SiC sheet, 0^2 ('^), 
which is derived from the calculated dielectric suscepti- 
bility, X = (^~l)/47r, multiplied by the distance be- 
tween two neighboring SiC sheets. ^41| Moreover, for the 
one-dimensional (ID) nanotubes, the measured optical 
spectrum is determined by the imaginary part of the op- 
tical polarizability per single tube, 0^2 (^)- This 0^2 is 
equal to the ef^^^- multiplied by the cross-sectional area 
of the supercell perpendicular to the tubular axis (A^^), 
i.e., a2 = (e§-^-) /4^. 



ELECTRONIC AND OPTICAL PROPERTIES OF 
SIC SHEET 

In the recent DFT-LDA calculations, a stable non- 
buckled graphene-like monolayer has been found for the 
isolated SiC layer 17, 42|. Despite the structural similar- 
ity, the unique massless Dirac fermion feature in graphene 
is no longer present in the SiC sheet because of the dif- 
ferent ionicities of the Si and C atoms. This heteropolar 
character induces an energy gap of 2.56 eV (the LDA 
value) at the K point of the hexagonal Brillouin zone, as 
shown in Fig. 2. The quasiparticle corrections to the 2D 
SiC sheet state are obtained through the GW calculations 
described above. The calculated quasiparticle band gaps 
and energy differenes at some high symmetry k-points are 
summarized in Table I. It is clear from Table I that the 
GW corrections to the LDA eigenvalues are rather sig- 
nificant, being more than 1.0 eV. We note that the GW 
correction to the LDA band gap of bulk 2H SiC is around 
1.0 eV. 43, 4^ Strikingly, the GW corrections transform 
the indirect LDA band gap {K M) into a direct band 
gap {K K) for the SiC sheet. The nature of the direct 
band gap in the 2D SiC sheet indicate that the lowest en- 
ergy exciton in the SiC sheet is optically active (bright). 
Furthermore, based on the simple zone- folding approx- 
imation, we could expect no momentum-mismatch in- 
duced optically inactive (dark) excitons below the lowest 
energy bright exciton in the SiC-NTs [45]. This anticipa- 
tion i is indeed confirmed by our GW+BSE calculations 
for both the SiC sheet and armchair SiC-NTs presented 
below. However, our GW+BSE calculations also show 
that this is not the case for very small zigzag SiC-NTs, 
as will be discussed in Sec. V. 

Figure 2 shows the LDA band structure (left panel) 
and also GW corrections to the LDA eigenvalues (right 
panel) of the 2D SiC sheet. We note that the optical di- 
electric function of the 2D SiC sheet is highly anisotropic 
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FIG. 2. (Color online) The LDA band structure (left panel) 
and GW corrections to the LDA eigenvalues (right panel) of 
the graphitic SiC sheet. The tt and a bands (left panel) are 
denoted as red and blue lines, respectively. In the right panel, 
the GW corrections for the tt and a bands are represented 
by red circles and blue diamonds, respectively, whereas the 
GW corrections of the NFE states are represented by green 
squares. The LDA E^~^* gap is marked as an arrow in the 
left panel. The valence band maximum is aligned at eV. 



(see, e.g., Ref. [13) • In particular, there is only weak dy- 
namical screening for the electric field perpendicular to 
the sheet in the energy range below 6 eV, whereas the dy- 
namical screening for the in-plane electric field is strong 
in this energy range (see Fig. 1 in Ref. There- 
fore, many-body interaction effects in the tt band (red 
line in Fig. 2) associated with the Pz orbitals, which ex- 
tend into the vacuum region from the sheet, would be 
less screened. Consequently, the tt band would have a 
larger quasiparticle correction and hence have a larger 
difference between the quasiparticle energies and LDA 
eigenvalues {E^^ — E^^^) (red circles in Fig. 2). On the 
other hand, because the in-plane a bonds are mainly con- 
fined to the 2D SiC sheet, electric screening effects on the 
a bands (blue lines in Fig. 2) are more significant, and 
hence GW corrections are smaller (blue diamonds in Fig. 
2). Furthermore, due to weak many-body interaction in 
nearly free electron (NFE) states which mainly lie more 
than a few As above and below the SiC sheet, the GW 
corrections to these conduction bands are small, being 
less than 0.3 eV (green squares in Fig. 2). This is similar 
to the NFE states in graphene[47j. These sophisticated 
screening effects enforce the notion that the quasiparticle 
calculations rather than the simple scissor-operator ap- 
proximation are necessary for low-dimensional systems 
such as the 2D SiC sheet. 

The excitonic effects on the optical properties of the 
SiC sheet can be examined by comparing the GW+RPA 



results with the GW+BSE calculations. Because of the 
huge depolarization effect in the 2-D planar geometry 
for light polarization perpendicular to the plane, we will 
focus on the optical absorption spectrum for light polar- 
ization parallel to the plane. Figure 3 shows the optical 
spectra from both the GW-I-RPA and GW+BSE calcu- 
lations (a) as well as the E^~^^ exciton wavefunction (b) 
of the SiC sheet. In Fig. 3(a), the first prominent peak 
located at = 3.25 eV in the GW+BSE absorption 

spectrum comes from a bright exciton due to the exci- 
tation between the tt and tt* states at the K point [see 
Fig. 2 (left panel)]. This is a strongly bound exciton 
with a large binding energy of Ef'^^^ = 1.17 eV as mea- 
sured by the difference between the E^~^* and the on- 
set energy of the GW corrected electron-hole continuum 
(^^TT-TT* ^ 4 42 eV). We note that the theoretical exciton 
binding energy in bulk 2H SiC is only 0.1 eV.[44] This 
shows clearly that the reduced dimensionality of a SiC 
sheet strongly confine the quasiparticles, and this signifi- 
cantly enhance the overlap between the electron and hole 
wave functions and hence the electron-hole interaction. 
In addition to this quantum confinement, the presence 
of the vacuum region reduces the screening and hence 
provides extra contribution to the large excitonic effect 
in the SiC sheet. An analysis of f \^s{'^e:'^h)\'^dTe and 
J |<l>5'(re, r/i)p(ir/i shows that the hole is mainly resided 
in the C-sublattice whereas the excited electron is on 
the Si-sublattice. In order to elucidate the extent of the 
^TT-TT* g-^(>^|^Qj^ wavefunction, we plot the electron ampli- 
tude square with the hole position fixed at the position 
slightly above the C atom (see Fig. 3b). Clearly, the elec- 
tron orbital is of the Pz character, and distributes mainly 
on the nearest neighbor Si atoms around the hole to form 
a fairly localized exciton. 

In the optical absorption spectrum of the SiC sheet 
a2{uj) from the GW+BSE calculations, the low-energy 
features in the energy range from 3.0 to 5.0 eV (Fig. 3a) 
are dominated by the tt ^ tt* transitions near the K 
point at the zone edge (Fig. 2), whereas the pronounced 
absorption peak at 5.83 eV is mainly due to the a ^ a* 
transition at the zone center. As shown in Fig. 3(a), 
the electron-hole interaction not only triggers a red-shift 
of the onset optical transition energies but also modifies 
their relative absorption intensities. Furthermore, be- 
cause the GW corrections transform the indirect LDA 
band gap of the SiC sheet into the direct band gap, as 
mentioned above, no dark exciton was found in the opti- 
cal spectrum until the photon energy up to 4.2 eV which 
is above the lowest bright exciton (£^{^~^*). 



QUASIPARTICLE BAND STRUCTURE OF SIC 
NANOTUBES 



An intuitive approach to construct a nanotube with 
a specific chirality is simply to roll up a layer of the 
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FIG. 4. (Color online) The fundamental band gaps of some 
zigzag (n,0) and armchair (n,n) SiC nanotubes as a func- 
tion of the tube diameter from the GW (solid symbols) and 
LDA [Tt^ (open symbols) calculations. The band gap of the iso- 
lated SiC sheet obtained from the GW and LDA calculations 
is also displayed as the solid and dashed line, respectively. 
Dotted lines between the symbols are a guide to the eye only. 



FIG. 3. (Color online) (a) The optical polarizability per unit 
area from both GW+RPA (blue) and GW+BSE (red) calcu- 
lations and (b) the E^~^* exciton wavefunction of the SiC 
sheet. In (a), the theoretical spectra are broadened with a 
Gaussian smearing width of 0.15 eV. In (b), the hole position 
(red sphere) is fixed at the top of a C atom (yellow sphere) 
and the squared electron amplitude (n = |$s(re,r^ = 0)|^ 
in an arbitrary unit (a.u.)) is mainly distributed on the Si 
atoms (cyan spheres) next to the hole. 



graphitic sheet along a specific lattice vector into a 
nanocylinder. This geometrical connection between a 
single-wall nanotube and the isolated sheet has inspired 
the appHcation of the zone-folding method to construct 
the zeroth-order electronic band structure and phonon 
dispersions for single- wall nanotubes [Uj. However, the 
zone-folding scheme alone cannot address the significant 
orbital rehybridization due to the curvature effect and 
also the many-body interaction effects in the nanotubes 
mentioned above. To include these effects, we have calcu- 
lated the diameter-dependence of the fundamental quasi- 
particle band gap of some small-radius zigzag and arm- 
chair SiC-NTs, as compiled in Fig. 4. First of all, we note 
that all the SiC-NTs are semiconducters, and this may 
be attributed to the different potentials of the Si and C 
atoms, being similar to the case of the BN nanotubes [isj. 
Secondly, an asymptotic convergence of the calculated 
band gaps of the zigzag and armchair SiC-NTs as a func- 



tion of tube size towards that of the isolated SiC sheet is 
observed in both the LDA and GW calculations. 

Interestingly, in contrast to the reciprocal radius de- 
pendence of the band gap in CNTs, we find a peculiar 
band gap reduction with the decrease of tube size for 
both zigzag and armchair SiC-NTs. This band gap re- 
duction of SiC-NTs with the increasing tubular curvature 
can be attributed to the curvature-induced hybridization 
between the tt* and a* orbitals near the conduction band 
minimum (CBM) in the small radius SiC-NTs, similar 
to that found in the very small diameter CNTs.[46] To 
illustrate this point, we display the charge density con- 
tours of the CBM of the zigzag (n,0) and armchair (n,n) 
nanotubes at the ID Brillouin zone center (F point) and 
zone boundary (Z point), respectively, in Fig. 5. Fig- 
ure 5 shows that in the small zigzag (n,0) SiC-NTs such 
as (5,0) and (6,0) nanotubes, a strong mixture of the tt* 
and cr* states in the SiC sheet now exists inside the tube 
to form a ring-like charge distribution, whereas the fea- 
ture of the Pz state is floating above the Si and C atoms 
outside the SiC-NT even under considerable structural 
distortion. Therefore, a substantial band gap reduction 
upon increasing the tube curvature can be expected. In- 
deed, this is confirmed by our LDA and quasiparticle 
calculations (see circle symbols in Fig. 4). In contrast, 
a much smaller degree of the orbital hybridization in the 
armchair SiC-NTs is induced by the folding-curvature, 
as revealed by the cylindrically symmetric distribution 
of the Si Pz component of the CBM at the Brillouin zone 
boundary (Z point) in the (5,5) nanotube (Fig. 5). This 
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FIG. 6. The energy differences between the quasiparticle 
eigenvalues (^^^) and LDA eigenenergies [E^^^) of the arm- 
chair (3,3) (upper panel) and zigzag (5,0) (lower panel) SiC- 
NTs. The corresponding LDA band structures are shown in 
the left panels for reference. The dashed lines denote the 
valence band maximum at eV. 



FIG. 5. (Color online) The charge density contours [p in 
10~^A ^)of the states at the conduction band minimum 
(CBM) at the ID Brillouin zone boundary (Z point) of the 
(3,3) and (5,5) SiC-NTs and also at the zone center (F point) 
of the (5,0) and (6,0) SiC-NTs. The Si and C atoms are rep- 
resented by the cyan and yellow spheres, respectively. 



gives rise to a minor size dependence of the band gap 
(square symbols in Fig. 4) for the armchair SiC-NTs. 
Nevertheless, a spectacular curvature effect in terms of 
strongly enhanced tt* — <j* hybridization does occur in 
the ultrasmall armchair (3,3) SiC-NT (see Fig. 5). This 
unusual orbital rehybridization effects on the band gap 
have also been recognized before in other ionic group III- 
V nanotubes[14, 45]. 

The quasiparticle corrections to the Kohn-Sham eigen- 
values from the GW calculations for two ultrasmall (ra- 
dius being ~ 5 A) armchair (3,3) and zigzag (5,0) SiC- 
NTs are given in Fig. 6, together with their LDA band 
structures. In the (3,3) SiC-NT, the GW corrections 
widens the LDA indirect band gap (2.13 eV) to the di- 
rect one of 4.09 eV, and also modify the band dispersions 
signicantly as a result of the complicated energy depen- 
dence of the GW corrections (see open circles in Fig. 6). 
Meanwhile, the small direct band gap (0.19 eV) of the 
(5,0) SiC-NT predicted by the LDA is increased to 2.02 
eV in the GW calculations. Furthermore, the GW cor- 



rections depend rather sensitively on the electronic state 
character, as illustrated in Fig. 6 for the (5,0) SiC-NT 
(open squares). We note that the curvature- induced hy- 
bridization between the tt and cr states in these ultrasmall 
SiC-NTs is responsible for the mixture of the quasipati- 
cle corrections of the otherwise orthogonal tt and <j states 
near the top of the valence band in the SiC sheet. This 
hybridization would become weaker gradually as the tube 
radius increases, and eventually diminish in the large 
tube radius limit where the self-energy corrections would 
be the same as that of the SiC sheet (as shown Fig. 2). 
The NFE tubular states are found in both (3,3) and (5,0) 
tubes as well as the small-radius chiral SiC-NTs. How- 
ever, these NFE states lie in the rather high energy re- 
gion (more than 5 eV above the band gap region from 
the LDA calculations) and hence play a minor role in the 
optical absorption spectrum. Therefore, they are not in- 
cluded in the present GW and GW+BSE calculations. 
Nonetheless, the NFE states can be shifted downward 
toward the band gap regime by doping electrons in the 
nanotubesfiit and consequently may significantly affect 
the optical properties of the n-type SiC-NTs. 



EXCITONIC EFFECTS IN SIC NANOTUBES 

In quasi-one dimensional nanotube systems, the wave- 
functions of quasielectrons and quasiholes can overlap 
significantly because of quantum confinement, and there- 
fore, electron- hole interaction may be considerably en- 
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TABLE II. The binding energies (Eb) of the first bright exci- 
ton (/i) and lowest energy dark exciton (Ki) below h for the 
single-wall zigzag and armchair SiC nanotubes (NTs) consid- 
ered here. The exciton binding energy (Eb) is calculated as 
the energy difference between the onset energy (E) of the cor- 
responding interband transition continuum and the excitation 
energy (Q). The tube diameter (D) and quasiparticle energy 
gap (E^J^) are listed for comparison. Symbol (i) denotes an 
indirect band gap in the armchair SiC-NTs. Moreover, the 
results for the SiC sheet which can be considered as the lim- 
iting case of an infinite diameter SiC-NT, are also listed for 
comparison. All the energies are in the unit of eV. 
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hanced. At the same time, the screening among the 
quasiparticles can be systematically changed because of 
curvature and reduced-dimensionality effects. For exam- 
ple, the phenomenon of anti-screening was discovered re- 
cently in the CNTs.[49] Therefore, the absorption fea- 
tures associated with the exciton formation can become 
prominant in the optical spectra in SiC-NTs. To study 
the excitonic effects on the optical properties of SiC-NTs, 
we have performed GW+BSE calculations by taking into 
account the electron-hole interaction kernel for several 
SiC-NTs of different sizes and different chiralities. 

In Fig. 7, the quasiparticle band structure and the 
imaginary part of the optical polarizability for the ultra- 
small armchair (3,3) SiC-NT are presented. Because of 
the strong depolarization effects in nanotubes [50], only 
the optical response for light polarized along the tube 
axis is significant. Therefore, only this polarization is 
considered here. For comparison, the optical absorption 
spectrum obtained without the electron-hole interaction 
in the GW+RPA scheme is also displayed in Fig. 7 (the 
blue-dashed curve in Fig. 7b). The onset energy (4.86 
eV) of this independent-particle optical spectrum corre- 
sponds to the lowest-energy dipole-allowed vertical tran- 
sition {E^'^) from the valence to conduction bands near 
2/3 of the T — Z line in the Brillouin zone which corre- 
sponds to the K point in the hexagonal BZ of the SiC 
sheet (see Fig. 7a). 

When the interaction between electrons and holes is 
included, excitons with specific excitation energies (ft) 
are formed to give rise to the prominent photoabsorption 
peaks. The binding energy of the lowest-energy exciton 




Photon Energy (eV) 



FIG. 7. (Color online) (a) The quasiparticle band struc- 
ture and (b) the optical absorption spectra of the armchair 
(3,3) SiC-NT. In (b), the optical spectra from GW+RPA 
and GW+BSE calculations are denoted as the blue-dashed 
and red-solid curves, respectively. The theoretical absorp- 
tion spectra were broadened with a Gaussian of 0.15 eV. The 
lowest-energy bright exciton (/i) is formed by the mixture 
of four interband transitions as indicated by the red arrows 
in (a). In the inset in (b), the squared electron amplitude 
(n = |$s(re,r^ = 0)|^ in an arbitrary unit (a.u.)) of the Ii 
exciton on a cross-sectional tubular plane with the hole fixed 
at the position of the red sphere (see text) are displayed. Si 
and C atoms are represented by cyan and yellow spheres, re- 
spectively. 



is defined as the energy difference between the continuum 
onset and corresponding excitation energy. For the arm- 
chair (3,3) SiC-NT, the optical absorption spectra from 
the GW+BSE calculations is dominated by the discrete 
exciton peaks [see the red curves in Fig. 7(b)]. The first 
bright exciton (/i) appearing at 2.86 eV with a large 
binding energy (£;§) of - 2.0 eV is formed by a mix- 
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ture of four pair of interband transitions, as indicated by 
the arrows in Fig. 7(a). These four different interband 
transitions couple strongly to make the electron ampli- 
tude more localized and highly asymmetric around the 
hole position located outside the tube but near the C 
atom, as displayed in the inset in Fig. (b), thereby re- 
sulting in a strongly bound character of exciton Ii. At 
the same time, the large electron-hole overlap in the (3,3) 
SiC-NT lifts the double degeneracy of the Ep^* exci- 
ton in the 2D SiC sheet to create the bright Ii exciton 
and a dark exciton with a higher (by 0.02 eV) excitation 
energy. Similar excitonic effects are also found in the 
larger armchair (5,5) SiC-NT (see Table II). However, 
the lowest optically active exciton in the (5,5) SiC-NT 
has a smaller binding energy of E^^ = 1.42 eV, because 
the reduced electron-hole overlap due to the decreased 
curvature. Moreover, apart from the bright bound and 
resonant excitons shown as red peaks in Fig 7(b), there 
also exist many dark excitons among the bright ones. 
However, we find no dark exciton below the Ii exciton in 
the armchair SiC-NT. As mentioned before, this can be 
anticipated from the fact that the isolated SiC layer has 
a direct band gap at the i^-point. In this case, there is 
no momentum mismatch when the SiC sheet is rolled up 
to form an armchair SiC-NT, and hence the possibility 
of the formation of the optically inactive lowest-energy 
exciton diminishes (iij. 

Next, let us concentrate on the excitonic effects on the 
optical absorption in the small radius zigzag SiC-NTs. 
In order to scrutinize the chirality effects, we performed 
GW+RPA and GW+BSE calculations of the zigzag (5,0) 
SiC-NT with a tube radius comparable to that of the 
armchair (3,3) SiC-NT. The quasiparticle band structure 
and optical polarizability 0^2 of the (5,0) SiC-NT is dis- 
played in Fig. 8 (a) and (b), respectively. In Fig. 8 (b), 
label Ij denotes the jth optically allowed excited state. 
The inset in Fig. 8 (b) is a contour plot of the electron 
charge density distribution of the Ii exciton with respect 
to a fixed hole position. Like the (3,3) SiC-NT, the pho- 
toexcited spectrum of the (5,0) SiC-NT is characterized 
by the discrete exciton peaks. The first bright exciton Ii 
is the most prominent feature in the spectrum, and con- 
sists of the different interband transitions near the zone- 
center F-point (represented by the solid arrow). As a 
result, the charge density distribution is not cylindrically 
symmetric. Similar results are found in the (6,0) SiC-NT, 
and the corresponding excitation energies are listed in 
Table II. Interestingly, although the (5,0) and (3,3) SiC- 
NTs have a similar curvature, the relatively smaller bind- 
ing energy (0.3 eV) of the h exciton in the (5,0) SiC-NT 
indicates a rather delocalized exciton. This pronounced 
chirality effect can be attributed to the considerable or- 
bital rehybridization in the zigzag SiC-NTs emphaszied 
above. This orbital rehybridization induces a significant 
band-gap reduction giving rise to the increasing of the ef- 
fective screening, and hence reduces the binding energy. 
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FIG. 8. (Color online) (a) The quasiparticle band structure 
and (b) optical absorption spectrum (broadened with a Gaus- 
sian width of 0.15 eV) of the zigzag (5,0) SiC-NT. In (b), the 
first bright exciton, /i, is mainly related to the interband 
transition indicated as a red-solid arrow in (a) whereas the 
transition by the red-dashed arrow in (a) corresponds to the 
third (/s) bright exciton in (b). The inset in (b) shows the 
squared electron amplitude (n in an arbitrary unit) of Ii with 
the fixed hole position marked by a red sphere. The Si and 
C atoms are represented by the cyan and yellow spheres, re- 
spectively. 



Furthermore, according to the zone- folding scheme for 
nanotubes, the lowest energy optically active interband 
transition at the K point of the 2D sheet can be folded 
into the F point to form a bright exciton with the lowest 
excitation energy in zigzag nanotubes. However, in the 
very small-radius zigzag SiC-NTs such as the (5,0) nan- 
otube, the strong curvature-induced orbital rehybridiza- 
tion has brought a higher conduction band down to the 
GW band gap (0.0 ~ 4.0 eV) region [see Fig. 8 (a)]. 
Since the optical transitions from the top valence bands 
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to this conduction band are dipole-forbiden, some low en- 
ergy dark excitons including the Ki dark exciton, appear 
in the small-radius zigzag SiC-NTs. As shown in Table 
II, this low-energy Ki dark exciton is also rather local- 
ized with a large binding energy, and this could make the 
zigzag SiC-NTs candidates for the tunable optical devices 
for temperature or external field sensors. 

Finally, we explore the extent of the exciton wavefunc- 
tion localization in the different SiC-NTs by comparing 
the intensity of their electron-hole wavefunctions in real 
space. In Fig. 9, the isosurface plots of the electron 
density with the fixed hole position (the red cross in the 
figure) of the first bright exciton (/i) of the armchair 
SiC-NTs suggest a significant local distribution in the 
direction along the circumference of the tube. Mixing of 
excitations to different subband pairs in the small-radius 
SiC-NTs is responsible for this tightly bound character. 
In addition, a spatial localization of the exciton wave- 
function along the tubular direction is also revealed by 
integrating out the electron coordinates in the perpendic- 
ular plane (red curves in Fig. 9). In contrast with the car- 
bon nanotubes, the deviation from the one-dimensional 
nature of the low-energy exciton is mainly attributed to 
the strongly mixing of the interband transitions. Similar 
tightly bound behavior of the lowest-energy bright exci- 
ton has also been discovered in the BN-NTsL24j- In the 
zigzag SiC-NTs, we find the similar anisotropic bound 
nature of the first bright exciton. The slightly broadened 
distribution of the integrated wavefunction intensity (see 
the blue curves in Fig. 9) suggests a less bound charac- 
ter in the zigzag SiC-NTs. The distinct spiky features 
in the curve are located on the cation (Si atom) planes 
and they indicate a charge transfer from an anion (hole) 
to its neighboring cations by photoexcitation. On the 
other hand, the pronounced peak structures are not ob- 
served in the armchair SiC-NT because of the specific 
atomic arrangement dictated by its particular chirality. 
Furthermore, the asymmetric distribution of with 
respect to the fixed hole position, is a consequence of the 
different nearest-neighbor cation-anion distances caused 
by rolling up the ideal symmetric 2D SiC sheet. This 
curvature-induced spatial symmetry breaking is also re- 
sponsible for the special nodal-like feature near the hole 
position in the smallest armchair (3,3) SiC-NT. Nonethe- 
less, by comparing the distribution of the exciton wave- 
function of the (3,3) SiC-NT with the larger (5,5) nan- 
otube, we note that this symmetry breaking will decrease 
gradually as the curvature decreases or the tube diameter 
increases. 



CONCLUSIONS 

In summary, we have employed the state-of-the art 
many-body GW and GW+BSE approaches to study the 
quasiparticle band structure and optical properties, re- 
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FIG. 9. (Color online) Isosurface plots of the electron distri- 
bution |$|^ of the h exciton (see Fig. 7 and Fig. 8) with the 
fixed hole position (indicated by a red cross) in the armchair 
(3,3) and (5,5) as well as the zigzag (5,0) and (6,0) SiC-NTs. 
The corresponding integrated intensities by averaging out the 
coordinates perpendicular to the tube axis, are also shown. 
Here the hole position is set at 0. The Si and C atoms are 
represented by the cyan and yellow spheres, respectively. 



spectively, of the SiC sheet and related SiC nanotubes. 
First of all, we find a direct quasiparticle band gap for the 
isolated 2D SiC sheet. Our GW band structure calcu- 
lations show that rather complicated orbital-dependent 
self-energy corrections are needed in obtaining accurate 
quasiparticle properties for the isolated SiC layer. The 
profile of the optical spectra is modified dramatically 
when the electron-hole interaction is included. In partic- 
ular, a strongly bright bound exciton with a large bind- 
ing energy (1.17 eV) is found to dominate the optical 
spectrum, because of the enhanced overlap between the 
electron and hole orbitals due to reduced-dimensionality 
effects and also the existence of vacuum region which re- 
duces the effective screening in the SiC sheet. 

Secondly, the quasiparticle band structure of the small- 
radius armchair and zigzag SiC-NTs is systematically 
studied within the GW approximation. Our detailed 
analysis on the charge density distributions reveals that a 
curvate-induced orbital rehybridization plays a vital role 
on determining the band gap of the small-radius SiC-NTs 
with different chiralities at both the LDA and GW levels. 
In particular, the quasiparticle band gaps as a function 
of tube diameter behave very differently, depending on 
the chirality of the SiC-NTs concerned. 

Finally, the calculated photoexcited spectra are con- 
sisted of discrete exciton peaks, thereby indicating strong 
excitonic effects in both the armchair and zigzag SiC- 
NTs. The optical absorption spectra of the small radius 
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armchair and zigzag SiC-NTs are dominated by the first 
bright bound exciton Ii with a significant binding energy 
up to ^2.0 eV. The highly asymmetric charge distribu- 
tion of the exciton Ii in the (3,3) SiC-NT is found to 
be a consequence of the strong coupling of the four op- 
tically allowed inter-subband transitions. Interestingly, 
we also find a quasi-zero dimensional bound character of 
the first bright exciton in the SiC-NTs and curvature- 
induced symmetry breaking affects the shape and size of 
this bound exciton. Moreover, our GW and BSE calcula- 
tions demonstrate that the simple zone-folding approach 
fails in predicting the low energy exciton characters in 
very small zigzag SiC-NTs, mainly because of the strong 
curvature-induced orbital rehybridization in these nan- 
otubes. We believe that the large excitation energy of 
~ 3.0 eV of the first bright exciton, with no dark exci- 
ton below it, may make the small-radius armchair SiC- 
NTs candidates for optical devices working in the UV 
regime. In contrast, the numerous dark excitons below 
Ii in the zigzag SiC-NTs may lead to potential applica- 
tions in tunable optelectric devices ranging from infrared 
to UV frequencies by external perturbations. 
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